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INTRODUCTION 


A dlffusor is a passive fluid device whose primary function is to con- 
vert the kinetic energy of a fluid into static pressure. It is a key element 
in centrifugal compressors and in power systems using expansion devices, 
in centrifugal compressors energy is imparted to the fluid initiaily as kinetic 
energy. The acceierated fluid is then siowed down through diffusing passages 
to increase its static prassure. In power systems using expansion devices 
such as turbines and magnetohydrodynamic generator channels the available 
exit pressure is fixed by the background. By attaching a diffuser to the 
end of the expansion device, part of the kinetic energy of the working 
fluid leaving the device can be converted into static pressure thus reducing 
the pressure at the exit of the expansion device beiow the available back- 
ground pressure. Since the efficiency of such power systems is a function 
of the pressure ratio across the expansion device, a decrease in the exit 
pressure of the device significantly increases the efficiency for a given inlet 
pressure. 

Considerable amount of experimental work has been carried out to 
determine the performance of straight two-dimensional diffusers [1,2, 3, 4]. 
Based on the results of these investigations Runstadler-Oean [4], and 
Reneau-Johnston-Kline [5] have generated diffuser performance maps to 
aid the designer. Attempts to theoretically predict diffuser performance 
have so far met with limited success, in order to better understand the 
diffuser flow and to predict diffuser performance when the design conditions 
(in particular diffusers with very large inlet blockage and heat transfer to 




the walls as encountered in magnetohydrodynamIc power generation applica- 
tions) are outside the range of avaiiabie experimental data there is need 
to develop theoretical methods to predict diffuser performance. Reliable 
theoretical methods to predict diffuser performance wiil also aid in exploring 
various ideas put forward to increase diffuser efficiencies. 

in the present work we have calculated pressure recovery coefficients 
aiong with other flow parameters for the diffuser geometries and the inlet 
flow conditions experimentaiiy investigated in Refs. [1,2, 3, 4). By comparing 
the calculated and experimental results we have derived an aigorithm to 
theoretically predict diffuser performance. The calcuiations were carried 
out using the method of Greywaii [6]. This method performs fuii viscous 
calculations, that is it caicuiates flow in a channei without separating it into 
a boundary layer and a core. Key points of the computational methodology 
of special significance to the diffuser calculations presented in this study 
are given in Appendix B. 
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COMPARISON OF IDEAL. CALCULATED. AND MEASURED Cp 
We follow the standard definition of pressure recovery coefficient. Cp. 

^ _ Pe-PI 

where p« and p| are the diffuser exit and inlet static pressures, and Pgj is 
the inlet stagnation pressure. 

Consider an incompressible and inviscid flow through a diffuser. We 
assume that at the diffuser entrance the velocity distribution normal to the 
flow is uniform (no boundary layer). Then, since the flow is inviscid. the 
velocity distribution will remain uniform as the flow proceeds along the 
diffuser. For such a flow application of continuity and Bernoulli's equaltion 
yields the ideal pressure recovery coefficient. Cp|. 



where AR is the ratio of the exit to the inlet diffuser areas. Later on AR 
will be referred to as the geometric area ratio. 

In Fig. 1 we have graphed Cp| as a function of AR. Also shown in this 
figure are Cpt< the pressure recovery coefficient we calculated, and Cp„,, the 
measured pressure recovery coefficient from Ref. [2]. Geometry of the diffuser 
and the inlet flow conditions used in the calculation of Cp( were identical to 
those used for the experimental Cp^. The flow conditions for Cpt and Cp„, 
are such that the flow can be considered incompressible. The diffuser length. L. 


Is equal to 12 times 2 Hq« where Hq Is half the distance at the inlet between 
the diverging diffuser walls. The boundary layer displacement thickness at 
the diffuser inlet, 6|, was equal to 0.04S times Hq. 

The difference between Cpi and Cpt is due to primary and secondary 
effects of friction. In the first place shear stress at the walls causes a 
loss in pressure. By summing up the calculated shear stress along the wail 
we have determined the loss in Cp| due to friction at the walls. This is 
shown in Fig. 1 by the broken curve marked Cf which is Cp| minus the Cp 
loss due to wall friction. The difference between the curve and the <=Pt 

curve is due to the fact that the area ratio seen by the flow is smaller than 
the geometric area ratio. This decrease in the effective area ratio seen by 
the flow is due to the blockage produced by the growth in nonuniformity of 
the veiocity distribution. 

As the flow proceeds down the diffuser the boundary layer grows 
causing larger and larger deviations from the uniform velocity profile 
assumed in the calculation of Cp{. The nonuniformity in the veiocity profile 
caused by friction is amplified by the rising pressure in the decelerating 
flow. That is the fast moving fluid along the outer edge of the boundary 
layer has to slow down less than the slower moving fluid near the wail to 
overcome the same rise in pressure. We shall call this, pressure rise- 
velocity profile-stretching. To exhibit the significance of this effect as 
compared to the effect due to friction, we consider the following comparison. 
The blockage calculated from theory at the diffuser exit, for the various 
area ratios considered, is shown in Fig. 1. The blockage, B, is defined as 
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wh«r« Ag is the geometric diffuser area and the effective aree seen 
by the flow. A^ff is caicuiated by dividing the mass flow rate by the product 
of the centerline density and velocity. We note, for later use, that for a 
channel with width much larger than its height, H, the blockage can be 
approximately written as. 



where 5^ is the displacement thickness. The exit blockage shown in Fig. 1 
is much larger than that caused by the frictional boundary layer growth 
alone (for fully developed turbulent flow between parallel plates the blockage 
is 0.125 if we assume the velocity varies as 1/7 power law). Thus the increase 
in blockage due to pressure rise-velocity profile-stretching is significant, 
it is also this increase in blockage caused by pressure rise-velocity profile- 
stretching that is responsible for the ever increasing difference bevween 
Cpi and Cpt as the area ratio increases, even though the loss in Cp due to 
friction at the wails decreases with increasing area ratio as shown in Fig. 1 
by the decreasing difference between Cpi and Cpf. 

The difference between Cp^ and Cp^ is due to partial flow separation 
that occurs in the real flow but is not accounted for in the theoretical calcula- 
tions. It is this partial separation that at first graduaily increases with 
AR and then eventually grows into a two-dimensionai stall that makes it diffi- 
cult to theoretically predict the diffuser performance. The separation pheno- 
menon in straight wall diffusers is different from that in flow past an airfoil 
at large angles of attack where the flow separates abrupty across the whole 
span of the airfoil. In straight wail diffusers (except for very large area 
ratio diffusers that have little practical interest) there is no abrupt flow 


saparation across the antira diffusar wall. Instaad ona talks In tarms of 
tha dagraa of flow saparatlon basad upon tha fraction of diffusar wall araa 
from which tha flow has saparatad. Culdad by tha axparlmantal data, Moora 
and Kllna {7], and Fox and Kllna [8] In thair pionaaring works davalopad 
tha concapts of various diffusar flow ragimas basad on tha dagraa of flow 
saparatlon. In our work wa hava adoptad thair classification. In particular 
the flow regime lines for no appreciable stall and large transitory stall where- 
aver indicated in our figures were taken from Fig. 2 of Raf. [5]. 

From tha proceeding discussion wa sea that if one is to theoretically 
predict the performance of a diffuser one must in the first place be able to 
determine from tha calculated flow properties whan a desired flow regime 
(no appreciable stall, maximum Cp, etc.) has bean reached, i.a. , a correia- 
tion between soma calculated flow parameter and the desired flow regime. 
And secondly, if possible, estimate the loss in calculated Cp( due to the 
partial separation. 

In the following discussion and in the figures, whenever comparing our 
calculations with tha experimental data reported in Refs. [1,2,3] wa have 
indicated blockage by its approximation 5^/H. This was done to conform 
to tha definition of the blockage used by these authors who carried out 
their diffusar experiments in large aspect ratio rectangular channels. 


CORRELATION OF DIFFUSER FLOW REGIMES 
WITH CALCULATED FLOW PARAMETERS 


Thtt first effort along theM lines was made by Ranaau [9]. Ha 
correlated the tine of no appreciable stall with d62/dx, the streamwise 
gradient of boundary layer momentum thickness. He found that for diffusers 
with small inlet blockage the line of no app/eciable stall is reached when 
d52/dx reaches a value of approximately 0.012. This criterion did not 
work for diffusers with appreciable inlet blockage. Our own computations 
confirmed Heneau's conclusions, in Fig. 2 we have presented Cp|, Cp(, etc., 
for diffuser flow with small inlet blockage; equal to 0.007. Values 

of experimental Cp„ were taken from Ref. [1]. Also plotted is our calcuia> 
tion of d62/dx at the diffuser exit. We find that dfi^/cix is indeed equal 
to about 0.012 at the line of no appreciable stall. However no such correla- 
tion was obtained for diffusers with large inlet blockage. 

All our efforts to correlate the boundary layer shape factor with flow 
separation failed. The separation criterion developed by Stratford [10] also 
did not yield correct results. The reason, we think, is the assumption in 
Stratford's analysis that close to the separation point the velocity profile 
very close to the wall varies as 


u 


Z 


.5 


where Z is the distance from the wail. Our computatons, as will be seen 
later, show that the partial separation in diffusers starts long before the 
exponent of Z reaches the value of 0.5. 


In our study ths point of maxifflum Cpg, is soisctsd for corrsiation. 

This is ths oniy flow rsgims that can ba unambiguously idantifisd from ths 
sxpsrimsntal data ovsr ths sntirs rangs of inist blockagss. Ths no appro- 
ciabis stall and largo tmasitory stall linos givsn in Rsf. [5] ars only fbr 
diffussrs with small inist blockags. For high inist blockags svsn ths lino 
for no spprociabis stall can bs past ths point of maximum Cp„,; sss for 
sxampis (an sxtrsms cass) graph for 6 i/Hq squai to 0.0SS in Fig. 8 whsrs 
ths lino of no apprsciabis stall is dsnotsd by a-a*. 

Ths corrslations prsssntsd in ths past havs boon bassd on ths global 
propsrtiss of ths boundary laysr. This was nscsssary sines ths boundary 
laysr calculations wars carrisd out using variations of ths momsntum intsgral 
approach. In ths prsssnt work sines full viscous calculations ars carrisd 
out ws havs sxpiorsd ths local (with rsspset to distance from ths wail) prop- 
srtiss of ths boundary laysr for a corrsla|ion paramstsr. 

To gain insight into ths dsvsiopment of ths boundary layer in diffussrs 
ws studied ths calculated velocity profiles at various distances along ths diffuser. 
In Fig. 3 is shown a typical set of such velocity profiles plotted in reduced 
boundary layer coordinates. 

u"*" = up/^T^ , and Z'*’ = , 

where is the shear stress at the wail, p the density, u the viscosity, 
u ths streamwise velocity, and Z ths distance from ths diffuser wail normal 
to the flow. In Fig. 3, 9 is the diffuser half angle and X the streamwise 
distance measured from the diffuser entrance. Dots on the curves represent 
the computational grid points along which u was computed. We note that 
the velocity profiles plotted on log-log scale ars not straight lines but have 


curvaturt. Tha curvatura incraaaaa with X. Similar curvas for othar 
valuas of d show that tha curvatura growa mora rapidly with X for largar 
valuat of 9. To bring out mora claarly tha curvatura of tha valodty 
profllas wa hava plottad in Fig. 4 tha gradiant 

d(in u’)/d(ln Z*) 

at a function of Z/H, whara H it tha haif diffutar haight at tha givan X. 

It can ba aatily tean that tha valocity dittribution 

u - Z" 

givet 

d(ln u*)l d(In Z"*") » n 

Thut wa can look at tha curvat in Fig. 4 at thowing variation of n with Z, 
if tha valocity dittribution is locaiiy (in terms of Z) approximatad by u«Z”. 

in our search for a parameter to correlate with maximum Cp we found 
that the value of n representative of tha vary near waii region velocity pro* 
file gave the best correlation. To be more precise the correlation parameter 
m presented in this work it the average value of n over the first four grid 
points (marked 2,3,4, and S in Fig. 3). This average was obtained from the 
slope of the tine drawn, using regression analysis, through the points 2,3,4, 
and 5 on In o'** versus InZ*^ curves. Calculations of m were carried out by 
putting additional statements in our computer program. 

For a given diffuser geometry and inlet conditions, we found that m 
increases monotonicaily with X, the distance from the diffuser entrance. 

The rate of increase of m with X increases with diffuser half angle and 
the inlet blockage, in Figs. S and 6 we have shown variations of the 
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paramtur m with X fbr diffarant valuaa of diffutar half angla. Figura S is 
for iniat biockaga, aqual to O.OOS and Fig. 6 is fbr <i/Hg aqual to 0.007. 

To datarmina tha corralation batwaan m and Cp^, maximum wa plottad m 
and Cpm varsus AR on tha sama graph. This was dona for diffarant sats 
of /aluas of L/2Hg and Typica* curvas of such comparison ara shown 

in Figi 7, 0 and 9. Also shown in thasa figuras is tha calcuiatad Cp| as a 
function of AR. Figura 7 shows tha comparison for diffusars of diffarant 
langths, L/2Hq, but all with tha sama iniat btockaga. i.a.. with 5 ^/Hq aquai 
to 0.007; tha axparimantai data for Cpm is from Raf. [1]. Figura 8 is for 
diffusars with diffarant iniat blockagas. but ail of tha sama langth. L/2Hg. 
aqual to 12; tha axparimantai data is from Raf. (2). Figura 9. onca again, 
shows tha comparison for diffusars of diffarant langths but ail with tha sama 
iniat blockagas aquai to 0. 12 (highast iniat blockaga invastigatad in this 
study); tha axparimantai data is from Raf. (9). 

From our studios of comparison of Cp||| and m varsus AR, wa conciuda 
that a valua of m aqual to 0.22 givas satisfactory corralation with Cpm maximum 
ovar tha antira ranga of L/2 Hq, araa ratio, and invastigatad. 

Using tha critarion that m aqual to 0.22 corrasponds to maximum Cp wa 
have thaoratically ganaratad locus of Cp maximum on AR varsus L/2Hg graphs. 
Rasuits ara shown in Figs. 10 and 11. Also shown in thasa figures ara pro- 
santly caiculr ‘ad and exparimantaily measured contours of constant Cp. 

Figura 10 is for tha smallest iniat blockaga, i.a., ^^/Hg aquai to 0.007, 
invastigatad in tiie prasant studies. Figura 11 is for tha highast iniat blockaga, 
aquai to 0.12, studied. In assessing how wall tha proposed criterion predicts 
tha araa ratio that corrasponds to maximum Cp for a given L/2Hg, piaasa note 
that in Figs. 10 and 11 tha point on a constant Cpm curve that corrasponds 


to CpiQ mtxifflum for a givon L/2Hp is tho point whora this curva is tangant 
to a constant L/2 Hq lino. 

For smaii iniat biockaga. tha iina of no appraciabia staii is waii aatabiishad 
(saa for axampia Pig. 2 of Raf. [S]). in our studios, wa aiso invastigatad 
corraiation batwaan tha vaiua of m and tha iina c'* no appraciabia staii. for 
iniat biockaga aquai to Q.007. Wa found that for this iniat biockaga vaiua 
of m aquai to O.iB corraiatod waii with tha iina of no appraciabia stall. 
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REDUCTION IN EFFECTIVE AREA DUE TO PARTIAL SEPARATION 

From our comparisons of Cp^ and Cpm versus AR curves, such as shown 
In Figs. 7, 8, and 9. we found that in the range in which Cpm it less than 
Cpm maximum, Cpm is well below Cpt at low inlet blockages, and that Cpm 
approaches Cpt as the inlet blockage increases. This implies that at low inlei 
blockage separation initiates well before Cpm maximum and moves towards 
Cpm maximum with increasing inlet blockage. Thus, at high inlet blockage 
one .i expect Cpt and Cpm to be nearly equal, in some of the cases 
with high inlet blockage Cpt was found to be smaller than Cpm< for Cpm 
less than Cpm maximum. Probably, the reason for Cp^ being less than Cpm 
is that the inlet boundary layer velocity distributions in the experimental 
studies were somewhat different than those modeled in our calculations. 

The flow in a diffuser is dependent on the inlet velocity distribution. Greater 
the boundary layer thickness at the inlet (i.e. , greater the inlet blockage) 
greater will be the impact on Cp of any differences between the experimental 
inlet velocity distributions and those modeled in our calculations. In our cal* 
culations the velocity distribution within the boundary layer at the diffuser 
Unlet was taken to vary as Z to the power 1/7. The boundary layer thickness 
was adjusted to meet the inlet blockage requirement. In the experimental 
setups the blockages was introduced by attaching lead channels of varying 
lengths upstream of the diffuser inlet to grow the boundary layer to give 
the desired inlet blockage. 

To provide a measure for the reduction in effective area (normal to the 
flow) caused by the partial flow separation we define. 
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A*' s ^ 

t yrrc 


and 


pt 


m 


4 ^ 


pm 


We note that Aj^ represents the effective area ratio seen by the computed 
flow« and A^ represents the effective area ratio seen by the actuai flow. 

Wo now introduce, 

A** aI' 

AR 

to represent the reduction in the effective flow area caused by the partial 
flow separation. 

In Fig. 12 we have plotted Z at the point of maximum Cpm as a function 
of inlet blockage, Bj. For inlet blockage greater than about 0.05 calculated 
Cp( is about the same as the actual Cpm at maximum Cpm> For inlet block- 
age smaller than 0.05, Fig. 12 can be used, along with the definitions of 
A^ and A^, to determine peak Cp^ from the calculated peak Cpf 

A systematic study of Z as a function of AR was ^rried out for the case 
of inlet blockage, 6 ^/Hq, equal to 0.007. This is the case most extensively 
studied experimentally [1]. The results of the calculation are shown in Fig. 13 
in the form of £ versus AR on log-log scale, for various values of L/2 Hq. Also 
shown in this figure is a line passing through area ratios corresponding to 
Cpm maximum. We note that irrespective of the diffuser length, £ grows smoothly 
as some power to AR. All the results shown in Fig. 13 are for flows below the 
flow regime "Fully Developed Two-Dimensional Stall" as given in Fig. 2, Ref. 

[5]. Results shown in Fig. 13 along with the definitions of A^ and a|^ can be 
used to determine Cpm that corresonds to a calculated Cp^. 

The remarkable observation shown in the results presented in Fig. 13 is 
that flow separaton in diffusers, at least below the flow regime of two-dimensional 
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Stall, is a gradual phanomanon avan through tha point of maximum Cp„,^ 

Thara is no abrupt changa in tha flow saparation pattarn associatad with 
Cpm maximum. Up to tha point of Cp„, maximum tha raduction in tha affactiva 
araa as saan by tha flow dua to tha flow saparation is lass than tha incraasa 
dua to tha incraasa in tha araa ratio and tharafora Cp^ incraasas. Bayond 
Cpm maximum a givan incraasa in tha araa ratio introducas a largar incraasa 
in tha flow saparatad area resulting in net decrease in the affactiva araa saan 
by tha flow and an associatad decrease in Cpm. 


CONCLUDING REMARKS 


Based on the results presented in the last two sections » one can construct 
various algorithms, depending upon the design specifications, to calculate the 
optimum diffuser geometries. For example, if the inlet conditions (inlet dif- 
fuser dimensions and the inlet flow conditions) and the diffuser length, L, 
are specified one can calculate the diffuser half angle (or the area ratio) for 
maximum pressure recovery as follows. Pick a value for the half angle, 6. 
and starting with the given inlet conditions calculate Cp( and m at the diffuser 
exit (i.e., at X equal to L). Repeat the calculations with larger and larger 
values of 6 until the value of m at the diffuser exit becomes equal to or greater 
than 0.22. Plot 9 and Cpt versus m. The value of 6 that corresponds to m 
equal to 0.22 is the desired half angle for the maximum pressure recovery. 

The value of Cp^ that corresponds to m equal to 0.22 is the maximum theoret- 
ical pressure recovery. From this Cp^ maximum, and the given inlet conditions 
we can now estimate Cp^ maximum using the results given in Fig. 12. 


APPENDIX A — NOMENCLATURE 
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L 

H 

W 

»0 

X 

z 

Ag 

Aeff 

Pe 

PI 

POI 

Cp 

Cpi 

Cpt 

Cpm 

AR 

a; 

A|. 

B 




diffuser length 

half the distance between the diverging diffuser walls 
half the distance between the parallel walls 
H at the inlet 

distance along the diffuser centerline 
distance from the diffuser wall normal to the flow 
geometric diffuser area 

diffuser area obtained by dividing the mass flow rate by the 
product of the centerline density and velocity 

diffuser exit pressure 

diffuser inlet pressure 

diffuser inlet stagnation pressure 

pressure recovery coefficient 


ideal Cp 
calculated Cp 
measured Cp 

diffuser geometric area ratio 

1/yi - Cpj , effective area ratio seen by the calculated flow 
1 // 1 - Cp^ , effective area ratio seen by the actual flow 


blockage * 1 - 


* TT 


for large aspect ratio diffusers 


streamwise velocity 

reduced boundary layer coordinates 


t 

I 





, i 

1 
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n a exponent of Z in u ** Z" 

m a average value of n in the very near wall region 
e a diffuser half angle 
6 q a boundary layer thickness 

6^ a boundary layer displacement thickness 

^2 - boundary layer momentum thickness 
a local wail shear stress 
p a density 

U = viscosity 
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APPENDIX B — COMPUTATIONAL METHODOLOGY 

Dlffus«r calculations ware carried out using the method of Ref [6]. 

In this method the flow in the channel is divided into finite streams and 
the finite difference equations are then obtained by applying conservation 
principles directly to the individual streams. The resulting finite differ* 
ence equations, solution technique, and the turbulence model used (based 
on the mixing length hypotehsis) are all given in Ref. [6]. Here we 
discuss two points of special significance to the present study. 

GRID SPACING 

Let be the wail. The next grid point Z 2 is seiected such that the 
iocai Reynold's number along Z 2 . i.e. 

kinematic viscosity 

is equai to 1200. This places Z 2 just outside the laminar sublayer and the 
transition zone. The grid point Zj is calculated from the relation 

Zj = Z2 + (Ho-ZjHj/N)^ 

where N is the total number of grid points (21 in the present calculations), 
and the exponent G is calculated from the requirement. 


EFFECT OF PARALLEL WALLS 


Lat W ba half tha distanca batwaan tha parallal walls. Tha affact of 
parallel walls was taken Into account by using instead of W in tha 
calculations. W^ is defined as 

Wg for 6g<W, and 

where 6g and 6^ are tha boundary layer thickness and the displacement 
thickness computed along tha diverging walls. 

For large aspect ratio (W»Hg), as for the Refs. 1, 2, and 3 cases, 
the proceeding correction to W Is not very crucial , although it was incorporated 
in our computations. However for an aspect ratio of one at the inlet, 
as for the Ref. 4 case, the proceeding correction is essential. 
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Figure 1 - Ideal « Cpj, calculated, St' and measured, Cp 

pressure recovery coefficients as functions of 
diffuser area ratio. Curve marked Cf is Cpj 
minus Cp toss due to friction at the wails. 

Also shown is the blockage, 6, at the diffuser 
exit as a function of AR. 
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Figure 7 ** Distributions of Cpt, Cp^, from Reference 1. 

and m as functions of AR for different diffuser 
lengths, L/2Hg. Inlet conditions: blockage = 

0.007; aspect ratio = 8; velocity = 69 m/s; 
fluid-air at N.T.P. 
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Figure 9 - Distributions of Cp^, from Reference 4, and 

m as functions of AR for different diffuser lengths. 
L/2 Hq* Inlet conditions: blockage = 0.12, aspect 
ratio = 1; Mach = 0.6, fluid-air. 
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Figure 10 - Measured and calculated contours of constant Cp for inlet 
blockage equal to 0.007. Also shown is the calculated 
locus of Cp maximum using the criteria m = 0.22 for 
Cn maximum. 
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Figure 11 - Measured and calculated contours of constant Cp for 

inlet blockage equal to 0.12. Also shown is the 
calculated locus of Cq maximum using the criterion 
m = 0.22. P 
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Figure. 13 - Reduction, due to partial flow separation, 
in the effective area seen by the flow, Z, 
as function of the diffuser area ratio, AR, 
for \ irious diffuser lengths. 





